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Abstract 

We are pursuing an experiment to measure the electric dipole moment of the electron using the 

O ■ 

^f^ . molecule PbO. This measurement requires the ability to prepare quantum states with orientation 






% 



of the molecular axis and, simultaneously, alignment of the electron spin perpendicular to this axis. 
It also requires efficient detection of the evolution of the spin alignment direction within such a 
state. We describe a series of experiments that have achieved these goals, and the features and 



C^ \ limitations of the techniques. We also discuss possible new approaches for improved efficiency in 

O \ this and similar systems. 
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I. INTRODUCTION 

A permanent electric dipole moment (EDM) of the electron, de, would violate both 
parity and time-reversal invariance [l|, since dg = 2deS (where S is the electron spin). 
There is substantial interest in measurements of dg with sensitivity beyond the current limit 



dp < 1.6 X 10 ^'' e-cm 



a. 



A non-zero value of de within the next few orders of magnitude 
would be a clear indication of physics beyond the Standard Model (SM). Moreover, a non- 
zero EDM within this range is predicted to occur in a wide range of theories that extend 



a 
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the SM |3|. 

Heavy polar molecules have a significant advantage over atoms in the search for d, 
a result of the strong hybridization of the atomic orbitals in such molecules, there is a large 
internal electric field £int = ^intn, where n is the direction of the molecule's internuclear 
axis. If de 7^ 0, unpaired electrons interact with this internal field giving rise to a linear 
Stark shift described by the Hamiltonian Hedm = ~^e ■ £-int- This leads to an observable 
energy shift in molecules that can be 2-3 orders of magnitude larger than what is achievable 
in atoms under typical laboratory conditions |5|, |6|] . 

Despite this advantage, molecules pose significant experimental challenges. For example, 
the thermal Boltzmann distribution limits the population in any particular rovibrational 
state. This reduces the counting rate in the experiment, and hence the overall sensitivity to 
dg. In addition, molecules suitable for measuring dg must have unpaired electron spins; such 
species are chemical free radicals and generally are not thermodynamically stable. This leads 
to substantial experimental difficulties and typically reduces counting rates even further. 



Measurements using the mo 



ecule PbO hold considerable promise for improved sensitivity 



to dg, as described in Refs. |7|, |8|, |9| and summarized here. PbO is not a free radical; rather, 
its ground state X(0)^S~'' has closed shells jlU]. The EDM measurement is conducted in the 
metastable a(l) pS+] state, which has two unpaired electron spins and hence is sensitive to 
dg. The a(l) state can be populated by laser excitation and has a lifetime of r^ = 82 yus. The 
i7-doublet substructure of this |i7| = 1 state allows it to be easily polarized with a modest 
external electric field £^ > 15 V/cm. Together these properties allow an EDM measurement 
using PbO contained in a closed cell, operating at substantial vapor density. These properties 
can in principle lead to an increase in the counting rate, over that available in molecular 
beam experiments, by as much as 7 orders of magnitude. The primary topic of this paper 



is the development of methods for achieving efficient state preparation and detection, which 
are needed in order to take full advantage of these unique properties of PbO. 

The f2-doublet structure of the a(l) state is central to the EDM measurement scheme 
using PbO, and introduces some new and unique aspects compared to previous EDM exper- 
iments. Here we summarize the salient features of this structure. Each a(l) vibrational state 
has a ladder of rotational levels with angular momentum J = 1, 2, ... and app roximate energy 



Ej = BrJ{J + 1), where the rotational constant Br = 271 x 7.054 GHz [10|, lll| for the t; = 5 
vibrational level used here. (We take h=l throughout.) Each rotational level contains a 
doublet of opposite parity states, separated by energy A^j = qJ{J + 1), where q <^ Br- The 
EDM measurement uses the J= 1 manifold of states, whose energy structure is shown in Fig. 
[T^. The J = l doublet is separated by Aq-^ = 27rxll.214 MHz [9\. These field-free states are 
designated by the ket \J, M, P, \Q\ = 1), where M is the projection of J along the laboratory 
2;-axis, P = ±1 is the parity, and Q is the projection of the electronic angular momentum Je 
along the molecular axis n: Q = J^. ■ n. Alternatively, these states can be written in a basis 
with signed values oi n, as \J,M,P,\n\ = l) = ^ {| J, M, fi = +l) - (-l)^P | J, M, fi = -l)} 



12l |. (Here and throughout the paper, kets without an electronic state label refer to the a(l) 
electronic state.) 

For the EDM measurements, external fields must be applied. An external electric field 
S = Sz leads to an ordinary Stark interaction described by the Hamiltonian Hst = —/J'a ' ^5 
here /x^ = Unh, where /la = 1.64 MHz/(V/cm) is the molecule- fixed electric dipole moment in 



the afl) state 



131] • Stark mixing of the f2-doublet states can be described by the Hamiltonian 



submatrix Hg^ of the two-state subspace (| J, M, P = — 1) , | J, M, P = +1)): 

2 J(J+1) / 

These states have their energy shifted above (higher state H) or below (lower state L) the 
field-free values; the magnitude of the shift 6st is given by 




The Stark mixing electrically polarizes the molecule, so that the resulting eigenstates 



have a nonzero expectation value of /i„. Defining the polarization direction as A^ = 
sign((^i^) ■ £), the resulting L (H) eigenstate has A^ = +1 (A^ = —1) for any finite 
value of £". In the limit of large £ (such that fiaS ^ '^fij), the f2-doublet states 
are fully mixed; to describe these states it is convenient to define the new basis set 
\J,M,N) = \J,M,n = N-sign{M)) = ^ {| J, M, P = +l) - (-1)^ | J, M, P = -l)}. These 
states have maximal polarization, {^i^) = Nfia\M\/[J{J + l)]i. Note that the M = sub- 
levels always remain unmixed, and that we have neglected mixing with much more distant 
states of different J. 

The effect of a magnetic field B = Bz on the system can be described by the effective 
Hamiltonian H^ = iXBdn li^)"^ ' ■^' where [xb is the Bohr magneton. This Hamiltonian 
leads to Zeeman shifts of sublevels with quantum number M by Sz^^ = Qjj i{£)iibBM 
compared to the field-free energies. The dependence of the (^-factors on £ is discussed 
in detail in section VI. Regardless of the value of £, roughly g\j' — Qj^ — jij\.i\ (here 
G\\ = 1.86 [9]), with differences at the level of a few parts per thousand. The relevant 
J = 1 level structure in the presence of magnetic and electric fields is shown in Fig. [H We 
assign special symbols to replace the subscripts H,L on g and 6z in some limiting cases: 
H{L) -^ P_{P^) when £ = (parity eigenstates), and H{L) -^ N^{N^) for large £ (fully 
polarized states). Since the paper focuses mainly on J= 1 states, we omit the superscript 
(J) for these states unless needed for clarity. Also, we sometimes refer to a generic Zeeman 
splitting as 6z, with no subscript. 

The desired EDM measurement sequence starts by populating a coherent superposition of 
the I J=l, M = ±1, A^ = +l) states and measuring the energy difference A7V+ between these 
two levels; or by populating the \J= 1, M = ±1, A^= — 1) states and measuring their difference 
A7v_ . Because of the different signs of {Smt) = £mtNz in these two states, the effect of de 
is of opposite sign in these two measurements, while the Zeeman shifts are nearly identical. 



(e//) , o^,., ,,„R U^^^ sri^ff) 



Explicitly, we write A^± = ±2de£>Z + '2gN±f^BB. Here £1^1^' ^ 26 GV/cm l]J,|l5| is the 



effective internal field acting on the EDMs of the pair of valence electrons in the a(l) state. If 
9n+ = 9n. exactly, the quantity Aedm = ^n+ — A^f^ = Ade£\nt i^ entirely independent of 
B. Hence, comparison of the energy shifts in the oppositely-polarized levels (with A^ = ±l) 
allows for the cancellation both of magnetic noise sources and of systematic effects due to 



B 



magnetic fields; i.e., the fi-doublet structure acts as an internal comagnetometer 

Suppose we initially prepare a particular initial state such as \ip{T = Q)) 
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FIG. 1: Energy level structure of the a(l) J = 1 manifold, (a) Structure with no applied electric 
field [£ = 0). The eigenstates have well defined parity. The f2-doublet splitting is Aq^ = 27rx 11.214 
MHz [9]. Zeeman shifts 5z due to the external magnetic field B are shown, (b) Structure with 
a fully polarizing electric field £ » An^Z/ia- Here sublevels with \M\ = 1 are states of mixed 
parity. Solid lines indicate the energy levels with only the Stark shift 5st taken into account. 
(Typical values are 8st ~ 27r x 40 — 100 MHz in our experiments.) Dashed lines indicate levels 
with both Stark and Zeeman shifts 5z included. (Typical values are 5^ ~ 27r x 50 — 500 kHz in 
our experiments.) Dotted lines indicate the additional effects of a non-zero dg. 

^(|J = 1,M = +1,A^ = + 1) + |J=l,M=-l,iV = + l)). After free evolution for 

time T, the state will evolve to |^(T)) = ■^{e-'^''^\J = l,M = +l,N = +l) + 
e*^/^ I J=l, M = — 1, A^ = +l)), where 6 = Ajy^T. Measurement of Atv^ is thus equivalent 
to measuring the time-dependent phase 9. Note that the state lipiT)) has, simultaneously, a 
well-defined molecular orientation (corresponding to the sign of A^) and a well-defined spin 
alignment perpendicular to the molecular axis. For the latter, we mean specifically that 
the expectation values (Jj) oc cos^ 6'/2 and ^J^) oc sin^^/2, while {J^ = {Jy) = (J^) = 0. 
(For precise definitions of orientation and alignment, see e.g. Ref. 16j.) In this sense, the 
evolution of \ip{T)) can be described as a precession of the spin alignment about the 2;-axis, 
and measurement of A7V+ can be achieved by monitoring the alignment direction. 

In order to implement the desired measurement, it has been necessary to devise new 
techniques to prepare this unusual type of state. This requires methods beyond simple laser 
excitation, since the molecule-oriented states (with definite value of A^) are not resolved in 
the presence of usual Doppler broadening. Moreover, these states are of mixed parity so 



that standard selection rules are modified. Similarly, the detection of A7v± also requires 
sub-Doppler methods, since under typical conditions in our work Aiv± ~ 50 — 1000 kHz 
(determined by the applied field B ~ 0.025 — 0.5 G). Our detection methods explicitly take 
advantage of the spin alignment precession, to directly measure 6. Although our discussion 
centers on the specific case of the electron EDM measurement using PbO, we believe our 
techniques for preparing states with simultaneous molecular orientation and spin alignment 
may find application in other systems and for other purposes (e.g. for encoding of qubits 



into the spin degree of freedom of polar molecules) |17 |. 

The discussion will begin with an overview of the key features of the apparatus used for 
these measurements. Next, simple schemes for state preparation and detection in the absence 
of an applied electric field, as used in our earlier work [9i|, will be described. In this context, 
the issues of population and detection efficiency relevant to the EDM measurement will 
be introduced. Then, three different schemes suitable for preparation of molecule-oriented 
states will be discussed. Finally, two additional detection schemes will be discussed along 
with their implications for improved sensitivity to dg. On several occasions we use these 
techniques to obtain various types of spectroscopic data on PbO; these measurements are 
described in the text as is relevant. 

II. APPARATUS 

The central apparatus for this experiment consists of a vapor cell consisting of an alumina 
rectangular frame, 2.5" tall with 3.5" wide sides (see Fig. [2]). The cell has 2" diam. holes 
covered by sapphire windows on four faces for optical access. The windows are held on the 



cell using gold foils as a sealant 18|. A stem tube, held at a temperature slightly below 
that of the main cell body, provides a reservoir of isotopically enriched ^°^PbO. A movable 
alumina plunger seals a hole in the top of the cell; this hole is used to evacuate the cell 
as needed. Electric fields are applied with gold circular electrodes (2" diam.) and annular 
guard rings (2.25" i.d., 2.5" o.d.), in a reentrant geometry to provide a suitable aspect 
ratio for achieving a uniform field (1.5" separation). Typically the guard ring voltage is 
^1.8 times larger than the electrode voltage to achieve optimum field homogeneity. The 
electrodes and guard rings are separated from each other and from the main cell by high- 
purity sapphire and/or beryllia insulating spacers. The cell is radiatively heated to 700° C 
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FIG. 2: (color online) Cross sectional top view of the apparatus. The apparatus consists of: (a) 
up to four layers of nearly cylindrical magnetic shielding; (b) the aluminum vacuum chamber; 
(c) heat shields; (d) the quartz and tantalum oven; (e) the alumina and sapphire vapor cell; (f) 
reentrant gold electrodes; (g) quartz lightpipes; (h) fluorescence collection optics and detectors; 
and (i) microwave electronics and horn. 



with resistive tantalum heaters supported by an opaque quartz frame. The cell temperature 
is measured with non-magnetic thermocouples, developed for this experiment, consisting of a 
Au/(50%Au-50%Pd) junction. The heaters, which are laid out to minimize their inductance, 
are driven with audio-frequency currents. The amplitude is modulated so that both drive 
currents and eddy currents are strongly suppressed during the measurement time after each 
laser pulse. The oven is surrounded by several layers of metal and quartz heat shields, 
and sits within an aluminum vacuum chamber. 2" diam. quartz rods are used as vacuum 
feedthroughs to convey light into and out of the oven to the cell. 

A pulsed Nd:YAG-pumped dye laser (rep. rate 100 Hz, energy ~ 15 — 25 mJ/pulse) is 
used to populate the a(l) \J = l) manifold of levels by driving the X(t> = l, J = 0) ^a(l)(f = 



5, J= 1) transition at wavelength A = 570 nm. The nominal laser propagation direction is 
defined as the y axis. We observe that the time-averaged laser spectrum can be described 
as the sum of two peaks with the same center frequency: one ~ 1.7 GHz FWHM containing 
~ 1/3 the power, and the other ~ 12 GHz FWHM containing the remainder of the power. 
The a(l) state molecules are detected by monitoring their decay fiuorescence along the 
a(l)(f = 5) — i>X(f = 0) channel at A = 549 nm. Fluorescence is collected along the ±a; 
direction by two of the quartz rods, which are polished on their circumference to act as 
light guides. The fiuorescent light passes through color glass and interference filters to 
block scattered laser light and blackbody radiation from the oven. The transmitted light is 
detected using photomultiplier tubes (PMTs) which are gated to further reduce scattered 
light. The signals are digitized, stored, and analyzed using a PC-based data acquisition 
system. 

III. PREPARATION AND DETECTION OF SPIN ALIGNMENT WITHOUT 
MOLECULAR ORIENTATION 

We begin with a brief review of the technique used to prepare and detect spin ahgnment 
in the |a, J = l) levels, when no electric field is present [.Qi]. The laser pulse driving the 
|X, J=1,P = +1) -^ |a, J=l) transition is linearly polarized along x. Standard electric 
dipole selection rules then ensure population of the state 

|V^(T = 0)) = |^i?l) = -L(|J=1,M = +1,P = -1) - |J=1,M = -1,P=-1)) (3) 

if the laser pulse duration ti satisfies the condition r^ <^ 5^^, so that the M = ±1 
sublevels cannot be resolved. This state evolves at time T to \ip{T)) = \iljp_(T)) = 
^(e-*'^^^|J = l,M = +l,P = -l) - e"^^^|J=l,M = -l,P = -l)). The measured fiuores- 
cence signal S(T) arises from decay of this state to all possible sublevels \J', M', P') of the 
X ^Ti~^{v = 0) state, where P' = (— I)'' for all levels. Then 

5(T)oc J]|(X,J',M',P'|e.r|^(T))|^ (4) 

/ 



here / denotes all possible final states of both the molecule and the polarization of the 
outgoing photon e. For detection along the x-axis, e = y or z spans the possibilities. 
Moreover, electric dipole selection rules ensure that only J' = and J' = 2 states contribute 
to the signal for |^(T)) = \iljp_{T)) (see Fig. [3]). Hence 

J' 

J'=Q,2M'=~J' j=y,z 

Within this sum, decays to M' = states are of most interest since here interference gives 
rise to quantum beats; hence we write S(T) = Sgh(T) + Sbg(T). Here Sgh [Sbg] is the quantum 
beat [background] part of the signal, given by 

Sgb{T)= J2 l(X,J',M' = 0,P' = +l|y|^p.(T))|' = 5j?[l + cos(25zT)] (6) 

J'=0,2 

and 

J' 

S'^9iT)= E El<X"^' = 2'^''^' = +l^.lV'p-(r))|' = S;;). (7) 

M'=-j' j=y,z 
AfVO 

Here we have used the fact that only e = y polarization contributes to SM = ±1 transitions 
that are detected along x. To derive the quantum beat signal it is useful to express y = 
— (r(+i)+r(_i))/-\/2 (where r(g) are spherical tensor components) and use the Wigner-Eckhart 
theorem to relate the matrix elements (X, J', M', P'| r^ |a, J, M, P) for different values of 
M,M'. We emphasize that 26z is exactly the desired energy splitting between the M = ±1 
components, and that the relative phase 6 = 26zT corresponds to the precession angle of 
the spin alignment in the x — y plane. 

We rewrite the signal as S{T) oc 1 + ccos{26zT), where c is the contrast of the quantum 
beats. The intrinsic value c = c^^^ = S^^ /{S^^^ + Sf^ ) is determined by the branching 
ratio for |a, J = 1,P=— 1) decays to |X, J = 2,P = +1) versus |X, J = 0,P = + 1) states. For 
allowed transitions, this could be calculated from Honl-London factors, but for the forbidden 
X[^S"*'] -^ a [^S+] transition no reliable calculation exists. Based on our measurements of 
relative line strengths for different J' -^ J transitions, we estimate c*^°^ ~ 50%. However, 
in practice there are other sources of background that reduce the observed contrast. In our 
experiments the broad spectral width of the laser leads to excitation of higher- J rotational 
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FIG. 3: Allowed decay channels for the a(l) J 
levels. Only final states with P' = +1 = (- 



1 states, a) Decays of the | J = l, M = ibl, P = — 1) 
are allowed, b) Decays of oriented molecule 
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states such as | J = 1, M = ±1, N = +1) . These mixed-parity states can decay to final states of both 
parities. In both cases, only decays to common (M' = 0) sublevels give rise to quantum beats. 
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FIG. 4: (color online) Typical Zeeman quantum beats, measured in the absence of an applied 
electric field {£ = 0). Both data (blue points, average over 70 laser pulses) and the fit function 
described in the text (red line) are shown. For this data, c = 9.2(1)%, uJb = 2tt x 299.8(1) kHz, and 
<j) = 5.02(1) radians. The background/signal ratio d/a = 1.10 is evident from the vertical scale. 

levels that contribute to the fluorescence signal but not the quantum beats. Under typical 
conditions, we flnd c ~10% as shown in Fig. HI 
We flt our data to a function of the form 



S{T) = aI{T) [l + c e-^/^" cos {ubT + 0)] + d. 



(8) 
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Here /(T) gives the temporal evolution of the signal in the absence of quantum beats, 
normalized such that /(O) = 1; typically /(T) ^ e~'^/'^^ (with Ti ^ 40 /is determined by a 
combination of the radiative lifetime r^ and coUisional quenching of the state on cell walls). 
The free parameters in the fit are the signal size a; the quantum beat contrast c, frequency cUb, 
and phase 0; the DC background d\ and the beat shortening time T(, (due to spin-decohering 
collisions). At various points in the experiment, I{T) was obtained by different methods such 
as a) deliberately applying an inhomogeneous magnetic field to rapidly dephase the beats; b) 
populating the |J = 1,M = 0,P = — 1) state by setting the excitation laser polarization along 
z] and c) filtering the quantum beat data using methods similar to those in Ref. [19|, l20|]. All 
methods give similar results and will not be distinguished from here on. The background 
d is due primarily to blackbody radiation from the oven. Typically the first ~ 10 yus of 
data is not fit: this data is contaminated by transients associated with PMT gating, and by 
fluorescence from excitation of the short-lived A(0) ^11 state by spectral wings of the laser. 

IV. STATISTICAL SENSITIVITY TO ENERGY SHIFTS 

Consider an EDM measurement using methods similar to those just described, but with an 
oriented molecular state prepared by some method. Then the beat frequency is u^, = Aiy±, 
and the EDM uncertainty 6de = Slui,/ i2£l^i '). Here ^cjf,, the uncertainty in the fitted 
value of ujf), determines the detectable size of d^. Hence, our goal is to minimize 6ujf) for 
measurements on oriented molecular states. A lower limit on 6ujf, is set by the shot noise 
on the fluorescence signal. The signal is S{T) = GN{T), where N{T) is the instantaneous 
photoelectron rate at the PMT cathode and G is an electronic gain. Straightforward analysis 
shows that the shot noise-limited uncertainty for a single laser pulse is 

5^W = ^ FjT^/n, d/a). (9) 

cT2y^iV(0)T2 

Here the beat coherence time T2 is given by Tg^^ = T^^ + T^^^ and the function 

F{Ti/Th,d/a) = J { i+T /t'' )^ "^ ^a- Under typical conditions T^ k, Ti; thus for a fixed 



molecular signal oc a, F ^ a/^ + Sq ^iid the frequency uncertainty degrades rapidly for 
background/signal ratio d/a > 1/3. The factor of \/2 in the numerator of Eq. [9] arises 
from treating the phase as a free parameter in the fit. (This analysis neglects the effect 



11 



of removing the first ~ 10/is of data from the fit.) Under typical conditions for the mea- 
surements discussed in sections III- VIII, A^ ~ 10^°/s and T2 ~ 30 fis. Over the course of 
these experiments d varied substantially, but typically d/a ?^ 2 — 10 and F ^ 4 — 9. With 
c ~ 10%, typically 600^ ^ 27r x 250 Hz under these conditions. In section X and Appendix 
A, the conditions determining A^ and d are discussed in detail. Our measured frequency 
uncertainty was typically within a factor of 1.0 — 1.3 of the shot noise limit given by Eq. [91 

V. QUANTUM BEATS IN FLUORESCENCE FROM ORIENTED MOLECULAR 
STATES 

In our initial experiments with oriented molecular states, we again used quantum beats 
measured in fluorescence to determine A^vi- In these cases, the only parameter that could 
be optimized to increase sensitivity was the beat contrast c. Our discussion in the next few 
sections thus concentrates on methods to improve c when detecting oriented states. In other 
experiments discussed later, we investigated alternate detection schemes in an attempt to 
further improve c and, at the same time, to increase both the signal/background ratio and 
the absolute signal size. 

Quantum beats in fluorescence from oriented molecular states such as {ipN+iT)) = 
('e-^^^+^/^|J = l,M = +l,Ar = +l)-e+^^^+^/2|J=l,M = -l,A^ = +l)VV2 are similar to 
those discussed earlier. However, the intrinsic contrast for these mixed-parity states is dif- 
ferent than for the P = —l state, since now decays to the states |X, J' = 1,M',P' = — 1) are 
also allowed. (This is the only molecular final state accessible from the P = +l components 
of |?/'7v+(T)).) A calculation similar to that above yields a slight reduction in expected in- 
trinsic contrast: c^^' ~ 45% for these polarized states. As before the observed constrast is 
substantially reduced due to fluorescence from other rotational states, but taking this into 
account our observations are in rough qualitative agreement with this prediction. 

VI. LASER EXCITATION WITH AN APPLIED ELECTRIC FIELD 

The simplest method to prepare states with simultaneous molecular orientation and spin 
alignment is essentially identical to the technique described above, but with the external 
electric field S present throughout the excitation and detection. A fully polarizing S causes 
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all four J = \M\ = 1 sublevels to be states of mixed parity, and the laser linewidth is 
significantly broader than the Stark shifted energy structure of the a(l) state; hence a x- 
polarized laser will excite a coherent superposition of all four of these levels. Initially, the 
excited state will be 

\MT=o)) = \^P_) 

= i(|J = l,M=-l,Ar = +l) - |J = l,M = +l,Ar = +l) 

_ |J=1,M=-1, A^=-l) + |J=1,M = +1, A^=-l)). (10) 

In ideal homogeneous fields S and B, the combined Stark and Zeeman shifts will cause this 
state to evolve as 



\MT)) 
- e"^^~"""'"--^'^|J = l,M=-l,iV = -l) + e"^^""""''"--^'^|J = l,M = +l,iV = -l)X.ll) 



Here us — Sst + Af7^/2, and we ignore the small effect of the electron EDM. Because 
1^5 ^ Szn± ^ V^2, even the small electric field inhomogeneity in the system {6S/S < 1%) 
leads to rapid decoherence between the N = +1 and A^ = — 1 components of \ip£{T)). The 
resulting system can be described as an incoherent mixture of two states, each of which has 
simultaneous molecular orientation and spin alignment: 

e"^'""+^|J=l,M = -l,iV = +l)-e+^'"^+^|J = l,M = +l,iV = +l) 

\vjn+[-I^ )) = y= ; (i/j 

e"''--'^\J=l,M = -l,N = -l)-e^'''--''\J=l,M = +l,N=-l) 
|^7v_(T)) = -= . (13) 

If gN+ = gN^, quantum beats in fluorescence from these states are indistinguishable. As 
described in our earlier work J9|, we initially expected this to be the case: the Stark-induced 
mixing of the J = 1,P = ±1 levels causes their gf-factors (initially different by the small 
amount gp^ — gp^ = 30(8) x 10^"^) to converge. However, our data indicates that the g- 
factors of the polarized states actually diverge as S increases (see Fig. |5]). This is explained 
below. 
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FIG. 5: Plot of difference in g-factors A^r between liigher and lower energy fi-doublet states, as a 
function of applied electric field E. Points are data; the dashed line is the result of the analytic 
calculation discussed in the text. The excellent agreement is evident. 

To obtain this data, we deliberately ran with a large magnetic field B {B ^ 0.5 G) to 
enhance the difference Ajv_|_ — A^v, — '^{gN+ ~gN-)f^BB. From the fitting procedure described 
above we extract the pure beat signal Sb{T), defined as the signal with the background d 
subr acted, the modulation envelope I(T) normalized away, and the beat offset subtracted: 
Sh{T) = ^jh. — 1. Under these conditions Sb{T) contains two Zeeman quantum beat signals 
of the same amplitude but nearly-equal frequencies; hence the term c cos {uj^T + (p) in the 
usual expression for S{T) [Eq. [8] must be replaced by c[cos (cj^+T + 0) + cos {uJb~T + 0)]/2 = 
ccos {oJhT + (f)) cos {AuJbT/2). Here uj^ = ("^6+ + ^6-)/2, and Aub = uj^^ — i^b-- The resulting 
pure beat signal is Sb(T) oc e~'^^^''[cos (uJbT + (p) cos {AuJbT/2)]. We determined AtUb = 
'^{gN+ — gN^)fJ'BB by finding the time Tq of the first node in the envelope function of Sb(T), 
such that cos (Au;fcTo/2) = 0. A signal of this type is shown in Fig. O 

The modification of the (^-factors of these states as a function of S is determined 
not only by mixing within the ^-doublet J = 1 levels, but also by mixing with the 
more distant J = 2 levels. We calculate guxi.^) as follows. We initially find eigen- 
states \iI)h,l{M)) of the two-state subspace (| J = 1, M, P = -1) , | J = 1, M, P = +1)) un- 
der the Hamiltonian H^^' (Eq. [1]). Next we include the effect of Stark-induced ad- 
mixtures of the J = 2 state. Since fia^ ^ B^ for all cases of interest, we use first- 
order perturbation theory to describe this mixing. The perturbed eigenstate is then 
\^h,l{M)) = \iJHAM))+V2+{M) \J = 2,M,n = +l)+r]2-{M) |J = 2,M,fi = -1), where 
r]2±{M) = {J' = 2, M', n' = ±1 Hst |V'h,l(M)) /(45^). Calculation of matrix elements of the 
form {J' = 2, M', Q'\ Hst |J = 1, M, Q) entering this expression is outlined in Appendix B. 
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FIG. 6: Quantum beats with an incoherent superposition of oriented molecular states. The data 
shown here have the background d subtracted and envelope function I{T) normalized away, i.e. 
they correspond to the pure beat signal Sh{T). The interference between the two beat frequencies, 
as described in the text, is evident. For this data, B = 540 mG, corresponding to the average 
beat frequency Of, = 1.35 MHz; the electric field amplitude £ = 118 V/cm; and we determine 
Aujb ^ 27r X 11.4 kHz. 



The first-order Zeeman shift 6zu ^ of this perturbed state, due to the magnetic Hamil- 



tonian Hmag, is given by 5zj, 



= {ipH,L{M)\Hmag\i^H,L{M)) and is calculated as follows. 
Off-diagonal (in J) matrix elements of Hmag are evaluated using the approximate form 
Hmag ~ G\\Je " T^B = G\\QB, whcrc G\\ = 1.86 is the component of the molecular G-tensor 
along the internuclear axis. Diagonal matrix elements of Hmag are calculated in the | J, M, P) 
basis, where (J, M, P\ Hmag \J, M, P) = gp±. This formulation is correct up to first order in 
both small factors Ag{0) and fiaS/Br. The (^-factor difference Ag{S) = gni^) — 5'l(^) is 
of primary interest, since it determines the zero-order effectiveness of the fi-doublet comag- 
netometer. The numerically calculated form of Ag{S) is shown in Fig. O where the good 
agreement with our measured values is evident. For the case fia^ ^ ^j^d a simple relation 
holds: 



Ag{S) = ^ Agio) + ^^ G^\. 



(14) 



flaS ^ ' 20 Br 

This method of simultaneously populating both states |V'Af+) and \ipN.) could in principle 
be used for measurement of d^. However, it requires working at large values of both B and £ 
in order to clearly resolve Acj^; this in turn increases the size of several systematic effects in 
the measurement of de- Moreover, we found it difficult to attain shot-noise limited sensitivity 
to Aujh from this type of data (although a more sophisticated data analysis method might 
remove this problem). For the moment, we have found it preferable to use a state preparation 
sequence that selectively populates either A^ = +lorA^ = — 1 states for a given laser pulse. 
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FIG. 7: (color online) Schematic of the laser-RF double resonance method for population of oriented 
molecular states. The laser pulse populates the | J = 1, M = 0, P = — 1) state. Then, a short RF pulse 
is tuned to drive the population to either A^ = — 1 (blue with frequency ujrfh) or A^ = +1 (red 
with frequency ujrfl ) • 

Two different double-resonance techniques were developed to achieve this goal, as described 
below. 



VII. ORIENTED STATE PREPARATION VIA LASER-RF DOUBLE RESO- 
NANCE 

The first double-resonance scheme for selectively populating one oriented molecular state 
worked as follows. Throughout the excitation and detection sequence S was present. A z- 
polarized laser pulse populated the \J = 1,M = 0,P = —1) eigenstate. Next, a short pulse of 
radiofrequency magnetic field /3^^ = Po cos {ujjiFT)y was applied to drive M = — > M' = ±l 
magnetic dipole (Ml) transitions to the oriented molecular states. When ujhf = uorf^ = 
6st {^rfl — ^st + ^Qi), the higher (lower) energy state \ipN_{0)) {\ipN+{0))) can be populated. 
Successful population of a \ipN±) state will result in appearance of quantum beats in the 
fluorescence signal. Note that in the limit of full mixing/polarization by S, the Ml transition 
strengths to the two oriented states are identical. The RF pulse durations must be sufficiently 
short to ensure that the linewidth of the pulses is broad enough to coherently populate both 
M = ±1 sublevels, as well as to cover the inhomogeneous broadening in the Stark-shifted 
average transition frequency due to inhomogeneity in the DC electric field S. This in turn 
requires that the RF field strength /3o be sufficiently strong to drive a vr-pulse within this 
short time. 

To apply the RF magnetic field, square single-turn Helmholtz coils were constructed 
around the cell body with their axis along the laser beam, in the y direction. The coils 
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were connected in parallel to reduce their inductance and driven via coaxial feedthroughs. 
Outside the vacuum chamber, it was necessary to impedance match the RF coils at both 
^rFh ^^'i ^rFl- This was achieved using two separate matching networks, one for each 
frequency, which were switched in or out of the line using relays. The coils were driven 
by a broadband pulse amplifier (1 kW maximum output power) excited with an arbitrary 
waveform generator for pulse frequency and envelope control. The RF field strength and 
homogeneity were measured using a calibrated miniature pickup coil, in a room-temperature 
mockup of the actual cell and heater geometry that reproduced all nearby conductive parts 
such as electrodes and heater elements. RF field strength as large as /?o ~ 1 — 1-5 G, 
corresponding to vr-pulse times T^p < 1.1 /is, could be applied at several different drive 
frequencies spaced by Af^^. The RF field was measured to be homogeneous to better than 
20% over the active volume of the cell, consistent with extensive numerical modeling of the 
system. 

Typical results for RF excitation are shown in Fig. [8l Here, uorf was fixed at either 
27r X 33.6 MHz or 27r x 44.8 MHz and the electric field strength £ was scanned so that 
\iPnJ) and \ipN+) came successively into resonance. Clear peaks in the quantum beat signal 
were observed at the expected positions. However, the maximum contrast achieved was only 
about 10-20% of that seen with direct laser state preparation as described previously. This 
poor RF excitation efficiency was inconsistent with expectations based on the measured RF 
field strength and homogeneity. Moreover, the width of the resonance peaks exhibited no 
sign of power broadening, which should have been marked at high RF power. 

In order to investigate this effect further, we studied the inverse process of driving M = 
±1-^M = transitions. Here, the states \i'N+) and H'n-) were initially populated using an 
x-polarized laser as in the previous section. In this case, driving the RF transition will result 
in disappearance of the quantum beats. (Note that for each of the two resonant conditions, 
each of the two initial states is resonant with a different final state | J=l, M = 0, P = ±l), 
so complete disappearance of the beat contrast should be possible.) As shown in Fig. [HI 
this effect was indeed observed, but showed substantially greater efficiency (approaching 
~ 70%) than in the beat appearance data (where the M = — * M' = ±1 transitions 
needed for preparation of oriented molecular states were driven). Moreover, the width of 
the disappearance resonances was consistent with the expected level of power broadening. 

The reason for the discrepancy between the apparent efficiency of excitation and de- 
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excitation in this system has not been fully understood. It appears to indicate that at the 
high drive fields expected to be necessary for a vr-pulse condition, the Zeeman coherence be- 
tween the \M\ = ±1 levels actually is being dephased more rapidly than the Rabi oscillation 
period. We suspect that this effect might arise from the mixed-parity nature of the \M\ = ±1 
states, as follows. The RF magnetic field (3j^p induces an RF electric field efip. Boundary 
conditions from the in-cell electrodes cause ejip \\ z over the active volume. This field has 
a node in the center of the cell, but at the edge of the electrodes (a radial distance d = 1" 
from the center) reaches the substantial value eRp{R = (i) ~ Pfip{dujfip/c) ~ 5 V/cm. This 
in turn induces a Stark modulation of the resonance energies ujrp^ ^^^ at the drive frequency 
ojRp. The depth of this modulation can be substantial, e.g. eRp{R = d)^a/'^ ~ 4 MHz, and 
is inhomogeneous at the 100% level over volume of the cell. It seems plausible that this field 
could lead to the observed dephasing, although no specific model we have applied appears 
to explain our results. 

As a result of the inefficient and poorly-understood results from RF preparation of the 
oriented molecular states, we investigated an alternate population scheme. This method, 
described below, uses microwave Raman transitions through the intermediate J = 2 rota- 
tional state. The much stronger electric dipole amplitudes associated with the rotational 
transition yield substantially greater excitation efficiency, and the microwave method has 
also proven easier to implement than the RF technique. 

VIII. ORIENTED STATE PREPARATION VIA LASER-MICROWAVE RAMAN 
DOUBLE RESONANCE 

As in the RF method, a z-polarized laser pulse populated the \J = 1, M = 0, P = — 1) 
eigenstate in the presence of S. Next, one resonant microwave pulse with polarization e^i 
and frequency cj^i excited these molecules to sublevels of the J = 2 manifold of states. 
Then, a second pulse with orthogonal polarization e^2 and frequency 07^2 = ^fii ~ ^st or 
^ti2 = ^iii + ^st + ^Qi stimulated deexcitation to the \ipN+) or \ipN-) state, respectively 
(see Fig. [9]). The linewidth F^^ of these microwave transitions is determined by Doppler 
broadening and/or inhomogeneity in 6st', under all conditions F^^ ^ ^ni,Sst so that the 
microwave pulses can easily resolve the states with well-defined molecular orientation. As in 
the RF scheme, the pulse durations must be sufficiently short to ensure coherent population 
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FIG. 8: Quantum beat signals with the laser-RF double resonance scheme. The quantum beat 
amplitude A (defined as A = ca) is plotted as a function of Stark shift 5st (a) Appearance of 
quantum beats following RF population of oriented molecular states. The black (red) trace is at 
high (low) RF power, as indicated in the legend. Even at the highest powers, the corresponding 
contrast is only about 10-20% of that seen with direct laser state preparation (corresponding 
to ^ ~ 3 on this plot). The lines broaden less than expected as the power is increased, (b) 
Disappearance of quantum beats following RF depletion of laser-excited states in the presence of 
£. Beat depletion as large as 70% was observed, and the width of these resonances was consistent 
with expectations due to RF power broadening. 

of the M = ±l states. 

We investigated these microwave Raman transitions both with e^i = x, e^2 = z (xz 
sequence); and with e^i = i, e^2 = £ (zx sequence). These sequences have different advan- 
tages and disadvantages. In the xz sequence, the first step creates an intermediate state 
\ipi{xz)) oc \J = 2, M = l) — |J = 2,M=— 1). After the second step, any residual population 
in this J = 2 superposition state can produce quantum beats that can complicate the data 
analysis (although these beats are at a different frequency than the desired J = 1 beats). 
Furthermore, Stark mixing in the \J = 2, \M\ = 1) states is incomplete over much of the useful 
range of S (both because A^.^ = 3Aq-^ , and because the Stark matrix elements are smaller 
between J = 2 sublevels). This means that the transition amplitudes in each step vary 
as £ changes, again complicating data analysis. The zx sequence (where the intermediate 
state is always \ipi{zx)) = \J = '^, M = 0, P = +l)) does not suffer from these difficulties. How- 
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ever, the second step of this sequence has a particularly small transition amplitude (due to 
unfavorable Clebsh-Gordan coefficients), making it difficult to saturate fully. 

Microwave generation starts with a two channel arbitrary waveform generator which 
outputs waveforms at a frequency of ~32 MHz. Each channel determines the excitation 
frequency for one step of the Raman process. At various points in our study of this process, 
we investigated vr-pulse excitation (by varying the area of the pulse envelope), as well as 
use of adiabatic passage (using linear frequency chirps). The ~32 MHz signal is multiplied 
by 32 times to ~1 GHz, filtered, and mixed with the 13.1 GHz signal from a phase locked 
dielectric resonator oscillator. A waveguide filter removes the low sideband leaving a signal 
at 14.1 GHz, which is doubled to match the splitting between the J = 1 and J = 2 states 
(nominally u^ ~ 27rx28.2 GHz) and amplified (up to 4 W output). A fast switch is used 
to direct the 14.1 GHz signal through either of two doubler/amplifier setups, which feed 
respectively into the x- or the z-polarized mode of an orthomode transducer whose output 
is connected to a conical horn. 

The horn output is located inside a 2" inner diam. Tefion tube, which slides over one 
of the 2" diam. quartz rods that extends into the vacuum chamber (ending 1 cm from the 
vapor cell). The Tefion tube and quartz rod act as a multimode waveguide, minimizing 
the angular spread of the microwave beam before it reaches the vapor cell. The microwave 
beam propagates opposite the laser beam, in the —y direction. The microwaves spread 
rapidly after exiting the light pipe, and undergo significant refiections from the sapphire cell 
windows; the power exiting the light pipe on the opposite side of the cell is reduced from the 
input power by ~ 10 dB. In addition, the amplitude of aj-polarized microwaves must vanish 
on the cell electrodes (parallel to the x — y plane) due to the usual boundary conditions. 
Hence, there is significant inhomogeneity in the microwave field amplitude over the active 
cell volume. Under most conditions the power delivered to the center of the cell was much 
greater than needed to saturate the transition; use of adiabatic passage rather than vr-pulses 
typically yielded somewhat higher excitation efficiency. 

The efficiency for driving microwave transitions was monitored by the appearance and/or 
disappearance of quantum beats in fiuorescence, as discussed in section VII. It was possible 
and often helpful to track quantum beats at the two different frequencies associated with the 
J = 1 and J = 2 states, by observing the power spectrum of the data. As mentioned above, 
the transition matrix elements varied considerably depending on the microwave polarization, 
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the value of £, and the specific J, M sublevels being addressed. Roughly speaking, however, 
at the maximum available power optimized vr-pulse transfer was achieved with pulse dura- 
tions of ~ 100 — 300 ns. Optimal adiabatic passage pulses were ~ 0.5 — 2 /is in duration, 
with a total frequency sweep of ~ 2 — 30 MHz (symmetric about the resonant frequency). 
Under a variety of conditions we have achieved > 50% transfer efficiency in single microwave 
pulses, and > 30% for the Raman transfer from \J = 1, M = 0, P = —1) to \ipNj,) or IV'at.)- 
Typical data is shown in Fig. [91 

We note in passing that the lineshape for (one-step, single-frequency) microwave excita- 
tion depends on the homogeneity of £ in the detection region, as a result of the different 
Stark shifts of the initial and final states. This effect made it possible to determine the 
r.m.s. variation 5£ /£ = 1%, in good agreement with expectations based on the distribution 
of £ found in finite element calculations with our electrode geometry. 

IX. ORIENTED STATE PREPARATION: LASER EXCITATION AND MI- 
CROWAVE ERASURE 

To date, the most efficient and robust method for creating states with simultaneous molec- 
ular orientation and spin alignment uses the technique we refer to as "microwave erasure" . 
This method begins as described in section VI: x-polarized laser excitation is used to pre- 
pare an equal- weight, incoherent mixture of \ipN+) and |^Ar_). Next, a microwave pulse, with 
polarization e^ = x, is applied for a long duration Terase ^ 5^^ so that the spin ahgnment 
precesses over many complete cycles during the pulse. The microwave frequency is resonant 
with either the \ipN+) ^ \J = 2,\M\ =2,N = +1) or the \iPnJ ^ \J = 2, |M| =2,N=-1) 
transition. The microwave power is chosen to be moderately high, so as to induce a tran- 
sition Rabi frequency Qr^^ such that 6z ^ ^r^ < Af^^. The combination of long pulse 
duration, large transition Rabi frequency, and microwave field inhomogeneity over the cell 
can lead to essentially complete dephasing of quantum beats from the J = 1 state that is 
in resonance with the pulse. (This "erasure" is similar to the behavior observed in the 
beat disappearance data for the RF transitions; however, here a single oriented state can 
be spectroscopically addressed as a result of the different f2-doublet and Stark splittings in 
J = 2 vs. J = 1 states.) The beats from the remaining oriented level thus have 50% of 
the maximum possible contrast, superior to any of the other methods tried. The complete 
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FIG. 9: (color online) Laser-microwave Raman double resonance state preparation. The plots 
show data from selective preparation of oriented states using microwave adiabatic passage. The 
r.m.s. amplitude spectrum [^^(wb)! of the pure beat signal S'fe(T) (see Fig. [6]) is plotted as a 
function of the beat frequency ujb- In both plots, the black solid traces show beat data with x 
laser polarization preparing an incoherent superposition of oriented states, as in Sec. VI. This 
provides a baseline for comparison to the case where oriented states are prepared. Preparation of 
oriented states began with z laser polarization, preparing the |J = 1,M = 0, P = — 1) state. Oriented 
states were populated with microwave excitation from this initial state, (a) Energy level diagram 
for the xz pulse sequence, (b) Results from the zx pulse sequence. The dashed blue trace shows 
the beat spectrum after both pulses, showing ~ 25% transfer efficiency. (No beats were observed 
after the first pulse.) For this data, the first (second) pulse had a frequency sweep centered 
around w^i = 27r x 28.220 GHz (a;^2 = 27r x 28.193 GHz), (c) Energy level diagram for the zx 
pulse sequence, (d) Results from the xz pulse sequence. The dashed blue (dot-dashed green) 
trace shows the spectrum after the first (second) pulse. After the first pulse, beats appear at 
the frequency corresponding to the Zeeman splitting between \J = 2,M = ±1) states. After the 
second pulse, ~ 45% transfer efficiency is achieved. The first (second) pulse was centered around 
w^i = 27r X 28.217 GHz, {lo^2 = 27r x 28.190 GHz). For both plots, £: = 39 V/cm. Each microwave 
pulse had 2 fis duration, during which the frequency was swept linearly over a range duj = 2tt x 11 
MHz. The second pulse arrived 0.5/is after the first ended. Both pulses had nominal power of ~ 1 
W. 

dephasing of the other state was verified in a variety of ways. For example, at high values 
of both S and B where the quantum beats of the two \4'n±) levels can be distinguished (as 
in section VI), the "erasure" of one beat frequency eliminates the interference between the 
beats. At lower values of S and B, the contrast is monitored as a function of microwave 
frequency and power. Near resonance and over a substantial range of microwave pulse power 
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FIG. 10: Confirmation of tlie microwave "erasure" method, (a) Plot of quantum beat amplitude 
A = ca versus microwave frequency, showing two dips at the expected resonant frequencies. The 
dashed horizontal line corresponds to the full beat amplitude seen with no microwaves applied. 
At each resonance, the beat contrast drops to approximately half its original size, indicating that 
the Zeeman coherence in one state has been eliminated. For these data, dB microwave power 
corresponds to ~ 1 W transmitted through the vapor cell, and the pulse duration was = 9/is. (b) 
Plot of beat amplitude A versus microwave attenuation, at the (resonant) microwave frequency 
indicated by the arrow in (a). Initially increasing the power increases the efficiency of "erasure" 
until the beats from one pair of states are completely erased and A has dropped by one half. 
Further power broadens the lines further, until the 2nd pair of states is also "erased" . Nearly-pure 
oriented states are prepared at the plateau indicating ~ 50% "erasure" in the range of powers 
-15^ -10 dB. 

and duration, the contrast after "erasure" remains at 50% of its value when the microwave 
pulse is not applied (see Fig. [10]). This indicates that a single oriented state has been 
selectively dephased. 



X. LIMITATIONS OF CURRENT METHOD FOR DETECTING SPIN ALIGN- 
MENT PRECESSION 

Next, we turn to alternate methods for detecting the precessing spin alignment in the 
oriented molecular states. The method of quantum beats in fluorescence has a number of 
features that induce serious limitations in the detection efficiency of useful signals, and in 
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addition give rise to large background signals. The effect of both factors on the ultimate 
uncertainty in measurement of d^ was noted in Section IV. In order to understand the 
motivation for the new techniques discussed here, it is useful to review the specific problems 
associated with the present technique. 

The overall efficiency of fluorescence detection is determined by several factors. The 
values of these factors varied over the course of the data described in this paper. Here we 
report the best (and most current) values: the solid angle of collection into the light pipes 
rjQ ^ 2 X 0.04 (with two light pipes); the efficiency of transmission through the optical 
filters rif ^ 0.1; the branching ratio for the detected vibrational decay channel rjER ~ 0.5; 
and the photomultiplier quantum efficiency rjge ^ 0.15. Both tjq and rif were determined 
by Monte-Carlo simulations and also checked by direct measurements. The total detection 
efficiency is thus r]tot = mVfVBRVg.e. ~ 5 x 10"^. 

The backgrounds in the current detection scheme are closely linked to the detection 
method and its efficiency. Since the a(l) state decays primarily to the X state, the detected 
fluorescence must be close in wavelength to that of the excitation laser. Detection at the 
laser wavelength leads to such large levels of scattered light that the PMTs (despite use of 
a proven gating circuit [21;]) cannot recover before the a(l) state has decayed. Hence we 
must use interference filters to select decays from the a(l) state at different wavelengths, 
corresponding to different vibrational final levels in the X state. Because the transmission of 
interference filters shifts towards the blue at non-normal incidence, scattered light remains 
a problem unless the selected fluorescence wavelength is to the blue of the laser. For this 
reason, we excite the X(f = 1) ^ a{v' = 5) transition at A = 570 nm and detect the 
a(f' = 5) -^ X(f = 0) transition at A = 548 nm using narrowband interference filters. 
(Note that at our cell temperature T ~ 700° C, the X{v = 1) state thermal population 
is roughly 1/3 that in the X(f = 0) state, thus reducing the excitation efficiency from its 
optimal value by the same factor.) The filter bandwidth (FWHM = 16 nm) is chosen as a 
compromise between wide-angle transmission of the fluorescence light (which improves with 
greater bandwidth) vs. increased levels of scattered light and blackbody radiation. Finally, 
the inability to resolve different rotational levels in the decay fluorescence leads to severely 
reduced beat contrast. 
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XI. SPIN ALIGNMENT DETECTION VIA OPTICAL DOUBLE RESONANCE 



These considerations led us to consider a new detection scheme. The basic idea, shown 
in Fig. [TT], is to probe the a(l) state spin ahgnment with a laser tuned to the a{v = 
5) — >C'(f' = 5) transition (at wavelength Aac' ~ 1114 nm) TJ, and subsequently detect the 
resulting C'(l) state population via its rapid decay fluorescence to the X state. This scheme 
has several advantages. The beat contrast c is expected to increase up to nearly 100%, 
since the double resonance method probes only the desired |a, J = 1,M = ±1) superposition 
state. In addition, the C'-X fluorescence occurs at blue wavelengths (see below), which leads 
to three improvements. First, the PMT quantum efliciency can be 3-4 times larger than 
at present; second, blackbody radiation is dramatically smaller; and third, scattered light 
from both excitation lasers should be trivial to eliminate, so that the optical fllter angular 
acceptance can be increased. The advantages are offset by a few difficulties. The most 
serious of these is related to the small cross section for the a— > C transition [7| : very large 
laser power (see below) would be needed to saturate the a— *C' transition over the vapor 
cell. In addition, the finite lifetime of the C state (re ~ 3.4 yus) acts as a low-pass filter to 
diminish the contrast of beats with frequency Ub > r^^/. Nevertheless, we project that the 
overall gain in sensitivity to de can be substantial using this detection method, and we thus 
have investigated it in some detail. 

Initial tests of this new detection scheme were conducted as follows. The state \ipp_) 
was populated as in section III, using an x-polarized laser with £ = 0. After a delay 
time T, the ensuing state \ipp-(T)) was probed by a second laser pulse, tuned to the Ql 
line ( J = 1 ^ J' = 1) of the a(l)[i; = 5] ^C'(l)[f' = 5] transition and with polarization 
^probe = £■ This results in population of the |C', J=l, M = 0, P = +l) state, with probability 
Pc' determined by the transition matrix element: 



Pc (X |(C',J=l,M = 0,P = +l|ep,„,e-r|V^P_(r)) 



|2 



e 



-iSzT 



(C',J=l,M = 0,P = +l|x|J=l,M = +l,P = -l) 



+e^'^^^(C',J=l,M = 0,P = +l|x|J=l,M = -l,P = -l))|V2 
oc l + cos{25zT). (15) 

Hence, the population of the C state refiects the direction of spin alignment in a(l) state 
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FIG. 11: (color online) Optical double resonance probe schemes using the a-C transition, (a) 
General features of the scheme, illustrated for the case of <? = 0. Here both excitation and probe 
lasers have x polarization, and a single sublevel of the C' state is populated, (b) Probe excitation 
of an oriented molecular state. Both P = +l and P = —l sublevels of the C' state are coherently 
populated by the probe laser, with their relative amplitude oscillating in time as the spin alignment 
of the |a, J = 1) state precesses. (c) Decays from the |C', J = 1,M = 0, P = ibl) superposition state 
created by probing oriented molecular states. Population oscillates between the P = +1 (dots) 
and P = —l (dashes) states. By collecting only l-polarized fluorescence (orange arrows), only the 
P = +1 population is detected and quantum beats with high contrast (c = 1) should be observed. 



(see Fig. [TT]). This C state population can be detected by monitoring its decay fluorescence, 
yielding a signal S(T) oc 1 + cos {26 zT) with 100% contrast (c = 1). 

For these measurements, the 1114 nm light was generated by a second Nd:YAG- pumped 
pulsed dye laser, downshifted to the infrared with a gaseous H2 Raman converter cell (Light 
Age). The pulses had duration of ~ 7 ns, bandwidth (averaged over many pulses) of ~ 5 
GHz, and power of ~ 0.7 mJ/pulse. To detect the C'-X fluorescence, a broad color glass 
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FIG. 12: Observed quantum beats using optical double resonance probe detection, with no applied 
electric field. The C'^X decay fluorescence signal is plotted as a function of time delay between the 
X polarized excitation laser [driving the X— >a(l) transition] and the x polarized probe laser [driving 
the a(l)— >C' transition]. As expected, near 100% contrast is observed. The DC background level in 
this data is believed due to light leakage through or around the low-quality blue-transmitting color 
glass filter used here. Subsequent measurements indicated that negligible blackbody radiation is 
present over the range of wavelengths corresponding to C'^X fluorescence. 



filter replaced the usual interference filter in front of a PMT, allowing for the detection of 
the decays from C'(t' = 5) to X(f' = — 9), corresponding to wavelengths Acx ~ 370 — 475 
nm. The variable time delay T was controlled with the computer used for data acquisition. 

Typical data from these experiments is shown in Fig. [121 Here, as expected the signal 
has (nearly) 100% contrast, validating the basic concept of this detection scheme. From 
measurements of absolute signal sizes we can deduce the strength of the a-C transition, and 
find it to be in reasonable agreement with estimates based on our previous measurement 
of the C -^a decay rate [Tj. This enables calculation of the saturation intensity for this 
transition (see below). We note in passing that spectroscopy of neighboring rotational lines 
in the a-C band allowed us to determine the C'{v = 5) rotational constant -8^=5 = 0.2395(3) 
cm^^ and vibrational energy G{v = 5) = 2625.1(1) cm^^ (taking T^ = 24947 cm~^ 10|])- We 
also determined an upper bound on the ^-doubling parameter qr/(v=F,) < 60 MHz. These 
are consistent with, but more accurate than previous determinations lOj. 

Next we discuss the application of this method to detection of states with simultaneous 
molecular orientation and spin alignment. Assume that such a state has been prepared, e.g. 
by any of the methods described earlier. Consider the signal resulting from probing the spin 
alignment of the state by excitation to the C state as just described. Because of the parity 
mixing in the oriented molecular state \4'n±), and because the f2-doublet splitting in the C 
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state is too small to resolve with the laser, the probe laser will excite transitions to both 
(unmixed, definite parity) M = states in C (see Fig. [TT]) . This results in population of 

the state iV'probe), where 

l^probe) OC ^ {C',J=l,M = 0,P\eprobe-1^\'4'P-iT)) 

p=±l 
OC e-'^^^ {C',J = l,M = 0,P = +l\x\J=l,M = +l,P=-l) 

+ e*^^^(C',J=l,M = 0,P = +l|x|J=l,M=-l,P = -l) 

+ e-^''^^(C',J=l,M = 0,P=-l|x|J=l,M = +l,P = +l) 

- e*^^^(C',J=l,M = 0,P = -l|x|J=l,M=-l,P = +l) 

OC cos((5zr)|C',J=l,M = 0,P = +l) 

+i sin {6zT) |C', J= 1, M = 0, P= -1) . (16) 

At first glance, it may appear that the information about the spin alignment preces- 
sion cannot be retrieved from this state, since the population in the C state Pc oc 
KV'probe I'j/'probe) P is independent of T. However, because the angular distribution of fiu- 
orescence from the P = +1 and P = — 1 states of C is different, quantum beats are indeed 
present in the fluorescence. A calculation of the beat contrast c similar to that in section III 
shows that in this case we expect c ~ 0.16 — 0.33, with the exact value again dependent on 
the unknown branching ratios for \C', J = 1,P = —1) -^ |X, J = 0, 2, P = +1) decays (see 
Fig. [TT]). A brief attempt to observe this type of signal was unsuccessful, but consistent 
with this expectation given the experimental conditions at the time. 

The quantum beat contrast for this type of signal (and all others discussed earlier) can 
be substantially improved by selecting a single polarization in the detected fluorescence. 
Suppose that only fluorescence with polarization eji = x{z) is detected. Then all AM = 
(AM = ±1) decays are eliminated from the sum in Eq. [S] and its analogues for the other 
cases discussed. In the present case of quantum beats following probe excitation, detection 
of only Sfi = z results in a signal with contrast c = 100%. This can be seen easily: the 
signal arising from population in the |C', J=l, M = 0, P = +l) state is out of phase with 
signal from the |C', J= 1, M = 0, P=-l) state, but the state |C', J=l, M = 0, P = +l) can 
decay only to states |X, J=l, M = ±l, P=— 1). (See Fig. [HI) Surprisingly, we have found 
(through both modeling and measurement) that the fluoresence polarization is preserved 
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with high fidehty even after transmission through the cyhndrical quartz hght pipes. Hence, 
it should be possible to obtain quantum beat signals with unit contrast using this method. 
In the current cell, the fluorescence polarization is scrambled by the birefringence of our 
sapphire windows, which leads to widely varying retardances over the angles of collected 
fluorescence. We are now inserting non-birefringent (YAG) windows and plan to look for 
polarization-analyzed signals in the near future. 

We envision using this detection method with a CW probe beam, since pulse-to-pulse 
fluctuations in the intensity of a pulsed probe beam would likely preclude shot noise limited 
measurements. Emerging technologies (diode laser oscillators and fiber and/or semiconduc- 
tor amplifiers) make it plausible to deliver sufficient power (P > 10 W) to saturate the 
a-C transition within r^. Under this assumption, we project a ~ 100-fold improvement in 
sensitivity to d^. Potential new sources of noise (e.g. due to laser intensity fluctuations) and 
new systematic effects (e.g. due to AC Stark shifts induced by the probe laser) appear to 
be controllable at the desired levels. We plan to implement this scheme in the near future. 

XII. MICROWAVE ABSORPTION PROBE OF SPIN ALIGNMENT 

Even under ideal conditions, the optical probe scheme described above suffers from the 
limited solid angle of fluorescence collection available (< 10% in our apparatus). In our 



a 



previous work (Ref. [7|), we pointed out that detection by monitoring absorption of a probe 
beam could in principle lead to much higher efficiency and sensitivity than any fluorescence- 
based scheme. The condition for this improvement is that the cell have column density of 
~ 1 for the probe transition. Unfortunately, there are no known electronic states of PbO 
which couple to the a(l) state strongly enough to meet this condition [7]. However, our 
success in microwave state preparation led us to consider using the J=l -^ J = 2 rotational 
transition (a fully allowed electric dipole transition) for absorptive probing of the J = 1 
state spin alignment. This can be achieved by monitoring the time-dependent absorption of 
x-polarized microwaves tuned to the \ipN±) -^ | J = 2, M = 0, P = +l) transition, in a manner 
exactly analogous to the optical probe described above. Note however that the final state 
of definite parity is easily resolved with the microwave probe because of its dramatically 
smaller Doppler broadening; hence there is no complication of the type discussed for optical 
probing of oriented states. 
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FIG. 13: Detection of |a, J = 1) state molecules by microwave absorption on the |a, J = 1) — ^ 
|a, J = 2) transition (with if^ = 0). The plot shows peak fractional absorption of the microwave 
power as a funtion of microwave frequency. Points are data (averaged over 512 laser pulses); the line 
is a fit to a Gaussian. For this data, the fitted amplitude is A^ = 3 x 10~^; the width is F = 27r x 52 
kHz FWHM, in reasonable agreement with the expected Doppler width Td = 27r x 42.5 kHz. 

In order to demonstrate the concept of this microwave probe method, we prepared 
an experiment to measure the a(l) state population by monitoring microwave "absorp- 
tion" in our apparatus. (In practice, we actually measured stimulated emission rather 
than absorption.) Unfavorable Clebsh-Gordan coefficients suppress the transition strength 
of the spin-alignment sensitive AM = ±1 transitions; hence, for this proof of con- 
cept we instead worked on the stronger AM = transitions by using ;z-polarized mi- 
crowaves. The states |J = 2,M, P = +l) were populated by laser excitation of the Rl 
line (|X, J = 1,P = —1) -^ |a, J = 2,P = +1)) with z polarization, then probed on the 
\J = 2,M,P = +1) —>-\J' = 1, M,P' = —1) transition. The apparatus was similar to that used 
in the microwave state preparation experiments, but here the transmitted microwave power 
was monitored as a function of time following the excitation laser pulse, using a crystal detec- 
tor (HP 422a). Electrical transients from the pulsed laser were subtracted off by periodically 
monitoring the signal with the microwaves far off resonance. 

As seen in Fig. [13], the peak microwave absorption is only A^ ^ 3 x 10"'* in our current 
configuration, even on this relatively strong probe transition. This is not sufficient for 
a competitive measurement of de using the microwave beam to probe the spin alignment 
precession. However, we believe that a few modifications to the apparatus could change this 
situation dramatically. For an absorptive probe measurement of this type, a long cylindrical 
vapor cell (rather than the roughly cubic cell now used) could increase the absorption column 
density. In addition, it should be possible to add a microwave Fabry-Perot cavity around 
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this longer cell, with partially reflecting mirrors designed to provide a finesse JF ^ 1, further 
increasing the effective column density. A detailed analysis of the signal-to-noise available 
with this technique is beyond the scope of this paper. However, we find that the primary 
noise source is likely to be thermal noise on the microwave detection amplifier; and that a 
factor of 1000 or more improvement in sensitivity to d^ could be obtained with this detection 
method, relative to the fiuorescence quantum beat method used at present. 

We note in passing that this microwave absorption probe method also enabled a variety 
of useful auxiliary measurements. For example, the transition amplitude and hence the 
absorption cross-section can be calculated accurately from the known molecular electric 
dipole moment /Iq. Hence, the absorption signal gives information on the absolute density of 
a(l)(J= 1) state molecules in the vapor cell, following laser excitation. Similar measurements 
on X state rotational transitions likewise yielded both the number of ground state molecules 
in the cell and the number excited into the a(l) state by the laser. These results are discussed 
in Appendix A. 

XIII. CONCLUSIONS 

In conclusion, we have developed several techniques for the preparation of molecular states 
with simultaneous molecular orientation and spin alignment, as well as several methods for 
detecting the precession of this spin alignment. These techniques were demonstrated with 
the a(l) state of PbO, as relevant for the measurement of de- Our experiments allowed 
a variety of auxiliary measurements, such as the dependence of the molecular (^-factor on 
the applied electric field £. The state preparation methods described here should enable a 
measurement of d^ with sensitivity within an order of magnitude of the current limit, and 
the improved methods for detecting spin alignment precession promise measurements with 
substantially improved sensitivity. 

We thank M. Kozlov for many useful discussions, including his derivation of the function 
g{S)\ and F. Bay, D. Kawall, and V. Prasad for their contributions to earlier parts of the 
experiment. We are grateful for the support of NSF Grant No. PHY0555462. 
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XIV. APPENDIX A: ABSOLUTE SIGNAL SIZES AND EXCITATION EFFI- 
CIENCY 

Here we briefly outline the methods for, and results of, determining absolute signal size 
and excitation efficiency in our experiments. The PbO X state vapor density is determined by 
microwave absorption measurements on the \X, J=l) —>■ \X, J = 2) rotational transition. We 
use the Stark-modulation method ^], with the modulation depth and transition strength 
determined by the known molecular electric dipole moment in the ground state, fix — 4.6 



D lOj. At our nominal operating temperature of 700° C, we find the PbO vapor density 
npbo (after accounting for the Boltzmann distribution over rotational and vibrational levels) 
to be npbo ~ 3 X 10^^ cm^^. We note that this is roughly one order of magnitude lower 
than the density assumed in our earlier work (Refs. [7|, |9|), based on the data in [23!]. A 
considerable literature exists on the partial pressures of various species such as PbO, Pb202, 



Pb404, etc. over solid PbO in this temperature range |24j. |25|. l26l. l27j Although there is some 
disagreement among various authors on the partial pressure of PbO, our result lies within 
the spread of data from these earlier measurements. 

We determine the efficiency of laser excitation from the X state to the a(l) state in 
three independent ways. In method 1, we use microwave absorption (as above) to mea- 
sure the absolute change in population of the \X, J = 1) state following laser excitation on 
the \X,J = 1) -^ \a, J = 2) transition (Rl line). In method 2, we measure the absolute 
population appearing in the \a, J = 1) state by the microwave absorption method discussed 
in Section XII. In method 3, we use the observed absolute size of our fluorescence signals, 
in combination with the calculated detection efficiency (see Section X), to determine the 
\a, J = 1) state population. All three methods are in agreement within a factor of 2, and 
indicate that the efficiency of laser excitation is rjexc ~ 2 — 4%. Note also that the agreement 
of method 3 with the others provides independent validation of our calculated detection 
efficiency rjtot = 5 x 10~^. 

The expected excitation efficiency can also be calculated as follows. As described in Ref. 
[3], the excitation cross-section aexc can be determined from the known a(l) state lifetime, 
the estimated Franck-Condon factors, and the laser and Doppler linewidths. (Note that the 
expressions for aexc in Ref. [71] contain a numerical error, and should be multiplied by a factor 
of 4 to obtain the correct value.) Based on this calculation, and assuming that ^ 1/3 of our 
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laser power is in a ?^ 2 GHz frequency band, we expect rjexc ~ 20%. However, this calculation 
predicts that rjexc should grow linearly with laser power up to our maximum available level. 
Instead, we observe instead that the \a, J = 1) state population grows roughly as the square 
root of the laser power, over an order of magnitude in power up to the maximum level. We 
believe this is because, while the spectrum of our laser is ~ 2 GHz wide when averaged over 
many pulses, the single-shot spectrum actually consists of 1-3 much narrower lines (each 
as narrow as ~ 30 MHz, the Fourier transform limit of our ~ 6 ns duration pulses). This 



behavior is known to be typical of pulsed dye lasers [2^ . In this case, for each laser pulse 
the small velocity classes of the Boltzmann distribution that in resonance with each narrow 
spectral line can be strongly saturated, while the other classes remain un-excited. Then, 
power broadening of the saturated signal leads to an increase in population proportional to 
the square root of the laser power. In this model, the expected excitation efficiency is reduced 
by a factor of 10-20 compared to the calculation assuming a single-shot linewidth of 2 GHz. 
This brings the calculated value of rjexc into better agreement with our measurements. 

XV. APPENDIX B: CALCULATION OF MOLECULAR MATRIX ELEMENTS 

For many of the calculations described in main text, it is necessary to calculate lab-frame 
matrix elements of operators defined in the body-fixed frame of the molecule. For exam- 
ple, the amplitude A for the |a, J=l, M, P = -l) -^ \a, J' = 2, M' = 2, P' = +1) microwave 
transition is determined by the off-diagonal matrix element of the molecular electric dipole 
moment operator fi^ = ^J^ 

A= (a, J' = 2,M',P' = + l|e^-^Ja, J=1,M,P = -1), (17) 

where e^ is the microwave polarization. Although there are a variety of texts which describe 



calculation of matrix elements of this type (see e.g. Refs. [12|, l29|), we believe the basic 
procedure may be unfamiliar to most atomic physicists, so we include it for completeness 
here. It is most convenient to work in the basis with signed values of Vt. Consider the matrix 
element M of a generalized irreducible tensor operator T^q (where k is the rank and q the 
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projection) between states in this basis: 

(J' k J 
-M' q M 

Here we have invoked the usual Wigner-Eckhart theorem: ( J'i7'| |Tfc| | Ji7) is a reduced matrix 
element and the factor in parentheses is a 3j symbol. This reduced matrix element can be 
expressed in terms of specific matrix elements evaluated in the molecule-fixed frame by the 
transformation 

im\\T,\\m) = V (-i)(-''-^V(2J' + i)(2J + i) I "^ ^ '^]{n'\n,'\n). (19) 

Here q' is a dummy index; at most, a single term with q' = Q — Q' contributes to the sum. 
The body-fixed states of the form 1^2) are defined in a frame (X, Y, Z) where the molecular 
axis n = Z, and Q = J^z- The matrix elements of the body-fixed operators used in this 
paper (n^ and Q) are trivial to evaluate in this basis. 
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